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HITTING TIMES OF RARE EVENTS IN BOUNDARY DRIVEN 
SYMMETRIC SIMPLE EXCLUSION PROCESSES 
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O 

Cn ' Abstract. In the boundary driven symmetric simple exclusion process con- 

}_^ sider an open set of density profiles which does not contain the stationary 

^ |- density profile. We prove that the first time the empirical measure visits the 

^^ • set converges to an exponential distribution. 

(N 

^ I 1. Introduction 

Ph . It has long been observed that in finite-state, reversible Markov processes the 

^^ I hitting time of a rare event is approximately exponentially distributed [IHl [131 [21 13] ■ 

For non-reversible dynamics much less is known. By estimating the total variation 
distance between the stationary measure and the quasi-stationary measure, Aldous 
[I] proved that the distribution of the hitting time of a rare event is close to an 
exponential random variable when the mixing time is small compared to the sta- 
tionary expectation of the hitting time. Fill and Lyzinski j23j proved that starting 
^ ' from the stationary distribution the hitting time of a configuration 77 can be repre- 

^— N . sented as an independent geometric sum of i.i.d. random variables if the probability 

QQ ' of hitting this configuration r/ at time t starting from ?/, viz. Pt(?7,?7), decreases in 

^^ . time. This representation permits to obtain bounds for the distance between the 

^b ' distribution of the hitting time and the distribution of an exponential random vari- 

f^ . able. Imbuzeiro ,25] proved that the hitting time of a rare event A is approximately 

CO ' exponential starting from a distribution 1/ if starting from v the probability of hit- 

ting A before the mixing time is small. Fernandez et al. [20] are presently working 
on this problem in the sequel of [6]. 

In this article we examine the hitting time of rare events in a well studied non- 
reversible dynamics, the boundary driven symmetric simple exclusion processes 
C^ . (BDSSEP). Beyond the complications arising from non-reversibility, this model 

presents a further difficulty in the lack of an explicit formula for the stationary 
measure. This obstacle is overcome by the use of a large deviations principle to es- 
timate the measure of sets, but prevents us from obtaining bounds for the stationary 
expectation of the hitting time with errors sharper than exponential. 

In the context of interacting particle systems the convergence of hitting times of 
rare events to exponential random variables has been abundantly investigated. Sev- 
eral results have been obtained for non-conservative dynamics, processes in which 
the local number of particles changes in time and which lose memory much faster 
than conservative ones. On the conservative side, which includes the the dynam- 
ics examined here, Ferrari et al. [3T] considered the case of a totally asymmetric 
one-dimensional zero-range process, and Ferrari et al. [2^ examined the case of the 
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one-dimensional symmetric simple exclusion process. This latter result was gen- 
eralized to any dimension and extended to independent random walks by Asselah 
andDaiPra glli]. 

The article is organized as follows. In the next section we state the main result. 
In Section |3] we present a general method to derive the asymptotic exponentiality 
of the hitting time of a rare event for finite-state, non-reversible continuous-time 
Markov processes starting from a measure not too far from the stationary measure 
in the sense of Lemmas l3.7[[3751 or l3.9] In Section 2] we estimate the expectation of 
the hitting time under the stationary state assuming a dynamical large deviations 
principle. In Section [5] we apply the results presented in the two previous section 
to the BDSSEP. 

2. Notation and Results 

The one-dimensional boundary driven symmetric simple exclusion pro- 
cess (BDSSEP). For A^ > 1, let Aat = {1, . . . , TV - 1}. Fix < a < /3 < 1 and 
consider the Markov process {ri^{t) : t > 0} on il]\[ = {0,1}'^" whose generator 



Ljv is given by 
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In this formula, rj — {ry(x), x G Ajv} is a configuration of the state space {0, 1}^" 
so that r]{x) = if and only if site x is vacant for ij; a'^''^r] is the configuration 
obtained from ry by interchanging the occupation variables ri{x), ri{y): 

and a^rj is the configuration obtained from rj by flipping the variable r]{x): 

ri{z) ii z ^ X , 



viz) 


a z^ x,y 


viy) 


if z = a; , 


T]{x) 


if z = y ; 



(cr''r/)(z) - . f \ -r 

1 I — ri[z) 11 z ^ X 

Hence, at rate a (resp. 1 — a) a particle is created (resp. removed) at the boundary 
site 1 if this site is vacant (resp. occupied). The same phenomenon occurs at the 
boundary x = N — 1 with /3 in place of a. 

Denote by -D(M+, S^n) the Skorohod space of paths from M+ to V,]y. Let P^, i] G 
Hn , be the distribution of the Markov process r]^ (t) when the initial configuration 
is 7]. The probability measure P^ is thus a measure on the path space D(IR+, fijv) 
endowed with the Skorohod topology. Expectation with respect to P^ is denoted 
byE^. 

The finite state Markov process f]^ (t) is irreducible and has therefore a unique 
stationary measure, denoted by i'^ p. The process is reversible if and only if a = /3, 
in which case the measure f^^ is a product measure. 
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The empirical measure. Denote by (•, •} the inner product in L2{[0, 1]) and 
set 

M:={peLoo([0,l]) : < p < 1} 

which we equip with the topology induced by the weak convergence of measures, 
namely a sequence {p" : n > 1} C M converges to p in M if and only if {p", G) -^ 
{p, G) for any continuous function G : [0, 1] -^ M. Note that M is a compact Polish 
space that we consider endowed with the corresponding Borel cr-algebra. 
Let d be a distance in M compatible with the weak topology, 

^(7,7') = 5]^|(7,i^fc)-(7',i^fc)|, (2.1) 

fc>i " 

where the continuous test functions F^ are absolutely bounded by 1. 

The empirical density of a configuration -q G ri^r, denoted by TT^{r]) £ M, is 
defined as 

where 1{^} stands for the indicator function of the set A. 

Denote by V the space derivative and by A the Laplacian. It has been proved 
in |16) that under the stationary state i'^ p the empirical measure tt^ converges in 
probability to the unique solution of the elliptic equation 

Ap = 0, 

p(0)=a, p(l)=/3. 



We denote the solution of this equation hy p — pa 



P- 



The dynamical rate function. To state the main result of this article we need 
to introduce the rate functions of the dynamical and the static large deviations 
principle of the empirical measure. We start with the dynamical one. 

For r > and positive integers TO,ri, we denote by C™'"([0,r] x [0,1]) the 
space of functions G: [0,T] x [0, 1] — > R with m derivatives in time, n derivatives 
in space which are continuous up to the boundary. We improperly denote by 
C"'"([0,T] X [0, 1]) the subset of C"'"([0,T] x [0, l]j of the functions which vanish 
at the endpoints of [0, 1], i.e. G e C™'"([0,r] x [0, 1]) belongs to C"'"([0,T] x [0, 1]) 
if and only if G(t, 0) = G(i, 1) = 0, t £ [0, T] . 

Let the energy Q : D{[0, T],M) ^ [0, oo] be given by 

Q{u) = 

sup I dt dxu(i,x)(VG)(i,x) - - / dt dxG(t,x)2 x("(i,x))| , 
G '-Jq Jq 2 Jo Jo -' 

where x ■ [Oj 1] ^^ ^+ is the mobility of the system, x(a) = a(l — a), and where 
the supremum is carried over all smooth functions G : [0,T] x (0,1) — >■ R with 
compact support. It has been shown in 10 that the energy Q is convex and lower 
semicontinuous. Moreover, if Q{u) is finite, u has a generalized space derivative, 
Vu, and 

1 r.. r.„(vu(i))^ 



Q{u) = 11 dt [ dii- 

^ Jo Jo 
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Fix a function 7 G M which corresponds to the initial profile. For each H in 

Co'^([0, T] X [0, f]), let Jh{-\i) ^ Jthi- D{[Q, T], M) — ^ R be the functional given 

by 

Jh{.u\^) ■■= {ut,Ht)~{7,Ho)- f dt{ut,dtHt) -\ j dt{ut,AHt) (2.2) 

Jo ^ Jo 

+ ^j dtVHtil) ~ ll dtyHtiO) - \j dt{xiut),{VHtf). 

Let /[o,T] ( • I7) : ^([0, T],M) — > [0, +00] be the functional defined by 
Ilo.T]iuh) := sup Jh{u\j) . 

HeC^^^([0,T]x[0A]) 

The dynamical rate functional /[o.t]('|7) • D{[0,T],'M) — !> [0,oo] is given by 

/roTl(^t7) = I ^["^^l^"'^) ifQH<oo , 
^ ' ^ [ cx) otherwise. 

The static rate functional. Denote by T^ : M — ;■ M_|_ the quasi-potential associ- 
ated to the dynamical rate functions /[o,t] • 

y (7) = mf inf {/[_T,o] {u\p) : u{-T) = p , u{0) - 7 } . (2.4) 

It has been proved in [U Theorems 2.2, 4.5 and A.l] that V is bounded, convex 
and lower-semicontinuous, and that V{p) > for all p ^ p. We are now in a 
position to state the main result of this article. 

Theorem 2.1. Fix an open subset 0/ M such that d{p,0) > and let Hq 
he the hitting time of the set 0, Hq = inf{t : tt^ {ri{t)) e 0}. Then, under v^ a, 
Hq/ EjjN [Ho] converges in distribution to a mean one exponential time. Moreover, 

if 

inf t/(7) = inlF(7), 

where represents the closure ofO, we have that 

hm i-logE,« [i/o] - inf F(7). 

7V->co TV °"5 7GO 

Finally, consider a subset 25 0/ M such that 

inf V{"f) < inf V{j) (2.5) 

where 23° stands for the interior of 23. Let Sjv = ('''^)^^('B) = {?? G ^n ■ 
Ti^ {rj) e 23} and let p^ be the probability measure on iljv defined by Pn{v) = ^{v ^ 
-^N^^a i3i'n)l^a H^^n)- Then, under p^i Hq/E^jN [Hq] converges in distribution 
to a mean one exponential time. 

This result holds in all dimensions, we restricted ourselves to dimension one for 
sake of simplicity. 
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3. Hitting times of rare events have exponential distributions 

Consider a sequence of irreducible, continuous-time Markov processes {tj (t) : 
t > 0}, N > 1, taking values on a finite state space J7jv- The points of Hn are 
represented by the Greek letters 77, ^. Denote by i^n the unique stationary state, by 
Lj\f the generator of the process, by Xj\f{ri), rj G SI at, the holding rates, by pn{ii, Oj 
S, ^ rj <E r^Tv, the jump probabilities, and by Rn{ii,£^) = Xn{v)pn{'>1iO the jump 
rates. In particular, for every function / : fljy — > M, 

Cef2iv 

We often omit the superscript N of 77^ (t) . 

For a subset A of 17 at, denote by Ha (resp. H^) the hitting (resp. return) time 
of a set A: 

Ha ■■= inf {s > : 77(5) € A} , 

Hj^ := inf {i > : ??(<) G A, 77(5) 7^ 77(0) for some < s < i} . 

When the set A is a singleton {77}, we denote i?{,,}, Hf , by H,j, f/'+, respectively. 

Let D(R_|_, Un) be the space of ^Ar-valued, right continuous paths with left limits 
endowed with the Skorohod topology. Denote by P^ = F^ , 77 G flw, the proba- 
bility measure on D(]R+, r^Ar) induced by the Markov process r]{t) starting from 77. 
Expectation with respect to P,, is represented by E^. For a probability measure fi 
in Qn, P^[ ■ ] = X^rjGOiv '^(^^ 'P''?! ' ]' with the same notation for expectations. 

Let Pt{ri,£^), t > 0, 7;, ^ G flN, be the semigroup associated to ri{t), Pt(?y,C) = 
P^[?7(t) = £]. Denote by \\fi — j^|Jtv the total variation distance between two proba- 
bility measures /i and v defined on SI at. Let T^"' be the mixing time of the process 

T^'^ = inf \t>0: max |1P,(77, • ) - z/w|1tv < 7} ■ 
Let An be a sequence of subsets of S^at such that 

lim miAw) = . (3.1) 

Denote by Hn — Ham the hitting time of ^at : 

Hm = inf{t > : 77^ (t) G An] , 
and by rAr(A^, An) the average rate at which the process jumps from A'j^ to An'- 






vn{A<j^) 



where Rn{S., An) = EceA„ -^^('^' 0- 

Nest statement is the main result of this section. It has to be compared with 
[TJ Theorem 1.4]. Instead of requiring that the mixing time is small compared 
to the the stationary expectation of the hitting time, we assume that the mixing 
time is small compared to the inverse of the averaged jump rate, rN{A%, An)~^ , a 
quantity easily estimated. Moreover, by |111 Lemma 2.3], for reversible dynamics, 
fN{A'j^, An)~^ is bounded by the expected value of the hitting time of An starting 
from the quasi-stationary state. 
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Theorem 3.1. Let An be a sequence of subsets offlN satisfying ()3.ip . Assume that 
T^"' <i; VNiAj^, An)^^ . Then, under vm, the sequence Hn /^v,^[Hn] converges in 
distribution to a mean one exponential random variable. 

Theorem 13.11 follows from Lemmas 13.41 and 13.51 We prove below in p.6p and 
Lemma 13.51 that 

\imudrN{A%,AN)^y^[H]s] > 0. 

In fact, under some assumptions this product converges to 1. To state this hy- 
potheses we need to introduce some notation. 

For two disjoint subsets A, B of the state space Q,n, denote by cap(A, _B) the 
capacity between A and B: 

cap(A5) = ^i^w(?7)Ajv(r/)P4ifB <i/i] . 

When the set A is a singleton, A = {77}, we write cap(77, B) for cap({?7}, B). 

Denote by {ri*{t) : i > 0} the stationary Markov process 7y(t) reversed in time. 
We shall refer to rj*{t) as the adjoint or the time reversed process. It is well known 
that 77* (i) is a Markov process on Q.^ whose generator L^ is the adjoint of Lm in 
L'^[vn)- The jump rates R^irj,^), 77 ^ ^ S ^n, of the adjoint process satisfy the 
balanced equations 

iyN{v)RN{v,0 = '^N{OR*Ni£.,v) ■ 

Denote by A*(77) = A(ry), 77 £ fl^, P*{'i]:0: V ¥" S, & i^N, the holding rates and the 
jump probabilities of the time reversed process ri*(t). 

As above, for each 77 € i^N, denote by P* the probability measure on the path 
space £'(M+ , 51 ^v) induced by the Markov process 77* (t) starting from 77. Expectation 
with respect to P* is denoted by E*. 

Lemma 3.2. Assume that there exists a sequence of subsets Bn, Bjq C A'j^, 

liniTv-i-oo '^7v(-B7v) = 1; such that 

hm sup Vl[HAr,<H„] =0, lim V 7^"^] ^ = «• (3-2) 

N^^^^B^ ^' " ^-'^^^^..j, cap 77, Aw 

Assume, furthermore, that 

limsuprJv(yl5^r, Ajv) < 00. (3.3) 

Af— >-oo 



The 



lim rAr(A^,AA,)E^„[ffjv] 



Proof. Fix 7; ^ An. By definition of the capacity, by equation (2.4) and Lemma 
2.3 in pij, and by the Markov property, 

cap(77,AAr) = cap*(r7,AA,) = ^ y^iO RN{iXWc[Hr, < HaA 

= Y. ^n{c)Rn[(:,o^i[h., <haA- 

€eA„ 

This sum is bounded above by VN{A'j^)rN{A'j^, An) < rjv(yl^, Ajv). On the other 
hand, if 77 belongs to Bn, by assumption p. 21) . the sum is bounded below by 
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{I — eN)i^N{A'^)rN{A'j^, An) > (1 — eAr)rAr(A^,ylAr), where eat is a sequence which 
vanishes as N '[ oo, and which may change from hne to hne. 
By [71 Proposition A.2], 

By the lower bound for the capacity obtained in the beginning of the proof and 
by (I3.2p . this expression is bounded above by 

In view of p.3p . this proves that 

hmSUprAr(A5^,ylAr)Ey„[i7Ar] < 1. 
AT-s-oo 

By p.2p and by the upper bound for the capacity obtained in the beginning of 
the proof, Ei,„ [Hn] is bounded below by 






»;Gi3jv 

This concludes the proof of the lemma. D 

Denote hy Nt, t >0, the number of jumps from A'j^ to An in the time interval 
[0,t]. Nt is a Poisson process and Mf, defined by 

Mt = Nt - RN{ri{s),AN)l{il{s) ^ AN}ds , 
Jo 

is a martingale. In particular, 

Note that {Hn <t} ^ {ri{0) G Ajv} U {Nt > 1}. Define 

Xt = 1{7?(0) e Aat} + Nt 

so that {Fat < t} = {Xi > 1}, and 

P.„[ifw < i] = F,„[Xt > 1] < E.jXt] 

< miAN) + tiyNiA%)rN{A%,AN) ■ 



(3.4) 



Lemma 3.3. Assume that T™"' -C rAr(yl^, Aat) -'^. Let jn, o'n be two sequences 
such that T™"^ ^ a^ ^ min{7jv, rjv(A^, Ajv)^^}. Then, for every t, s > 0, 



■ vn 



{Hn > (t + s)7Ar] - P,„ [i^AT > S7Ar] P.„ [i?Ar > t7iv] 
< 2j.Ar(AAr) + 2 ^at rAr(yl^^, Aat) + (1/2^"/^"". 
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Proof. In view of the definition of Xt, we have to estimate the difference 
Let (Tat be a sequence such that T^'^ ^ cat ^ rAr(A^, Aat)^^. Clearly, 



[XsfK — Oj - Pi,„ [Xs^^ - Xct„ — oj 



< 



[Xa. > 1] 



and, by (|3.4p . this last probability is bounded by i/Ar(AAr) + aN rN{A'^, An)- By 
stationarity, a similar bound holds for the absolute value of the difference 



It remains to estimate the absolute value of the difference 



= 0]. 



[Xp 



■ vn l^t7JV 



0, X 



{t + s)^N 



X: 



t^N+CN 



0] 



By the Markov property, this expression is equal to 

^{^t-fN = 0} |P,,(t^„) [Xs^^ - X„^ = Oj - Piyjv [^SfN " ^CTN = Oj I 

This expectation is absolutely bounded by 



sup 



E, 



^[^. 



'?("■«) l-^siN—crN 



OJ 



UN [^S7JV-CTiV 



Oj < (1/2) 



o-n/T" 



where we used the definition of the mixing time in the last inequality. This concludes 
the proof of the lemma. D 



Let 0N be given by 

On = inf |t > : P^„ [Hn > t] < e^^j 
Note that V^„ [Hn > On] < e-\ Hence, by (lO]) . 



(3.5) 



1 



< 



[Hn<9n] < i^NiAN) + 9NiyNiA'N)rNiA%,AN) 



Since j/7v(^Af) vanishes, we deduce from this inequality that 



\iTnmi9NrNiA'ij,AN) > 

N 



(3.6) 



In particular, 0n ^ T- 



mix 

N ■ 



Lemma 3.4. Assume that T™"' ^ rN{A'^^,AN)~^- Let On he the sequence de- 
fined by (|3.5p . Under vn, the sequence of random variables Hn/9n converges in 
distribution to a mean one exponential random variable. 

Proof Since T^'"" < rjv(A^, ^jv)"\ by dSJl), ?!"'' <t: On- By Lemma EJ with 
JN — On and some sequence aN, T^^'^ ^ aN ^ I'NiAj^, An)^^, we have that 

Hn 



lim 

AT-j-oc 



> t 



^N 



t>0 



D 



Lemma 3.5. The sequence On introduced in (13. 5p satisfies 

IEi/„ [ Hn ] 



lim 

A— i-oo 



1 . 



7A 
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Proof. Let 

ONiv) ■■= inf {t > : P,Ji/w > t] < e"^} , rjeflN , 
and let 9^ — max,jgj2„ 6j\j(ri). We first claim that 

lim 9n/0n = 1 . (3.7) 

It is clear that 9/^ < ^jv • Indeed, ii t > 9j^i, t > 9^ (»/) for all rj £ J7jv , so that 

¥,JHn >t] = J2 ''n{v)V^[Hn >t] < e^i . 
i)sr2jv 

Hence, 9^ <t and Sjv < On- 

To prove the converse inequality, let 6]^ (a), a > 0, be given by 

9N{a) := inf {i >0 : F^^[Hn > t] < e^'^} . 
For any ij G Hn, e > 0, L > 1, 

F^[HN>9Nil + e)+LT'j^'^] < E^[f,^^lt^^...-,[Hn > 9n{1 + e)] 
By definition of the mixing time and of 0Ar(l + e), the last expectation is bounded 

by 

provided 2-^ < e-^[l - e"^]. Hence, 6'jv(?7) < 6'Ar(l + e) + LT^J"'^ for ah 77 e 17 a, so 
that 9n < 6'jv(1 + e) + LT™'^. 

Denote by i?jv the right hand side of the inequality appearing in the statement 
of Lemma 13.31 with 7jv replaced by 9^. Iterating k — 1 times this estimate, we 
obtain that 



P,^[HN>9N/k] < (P,„[FAr>0jv] +fci?JV^ 

Applying once more Lemma 13.31 we get that 

F.^[HN>{k + l)9N/k] < P,^[Hn > 9n]P.^[Hn > 9N/k] + Rn , 

so that 

P^^[HN>{k + l)9N/k] < e-^{e-^ + kRNY'^ + Rn ■ 

Since Rn vanishes, if fc > e^^ this expression is bounded by e"^^"'"'^-' for A'' suffi- 
ciently large. Therefore, 9m{^ + e) < (1 + k~^)9N for all N large enough if A: > e~^ . 
Taking k = [e^^] + 1, where [a] stands for the integer part of a, we conclude from 
the previous two estimates that for N large enough 

^A' < {i + j-^)9m + LT^"- 

provided 2~^ < e^^[l — e~'^]. This proves that for every e > 0, lim sup^ (6* at 76* at) < 
1 + ([e^^] + 1)^^, i.e., that limsupAr(^Af/^Af) < 1, proving claim p. 71) . 

It follows from Lemma [3^ and p.7p that Hm/9m converges in distribution to a 
mean one exponential random variable. We claim that 

lim ^M^hd ^ 1 . (3.8) 
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To prove p.8|) . we change variables to obtain that 

^-^^ [Hn] = 6]^' P,AHN>t]dt= F,^[HN/ON>t]dt. 

Jo Jo 



/o 



It remains to obtain a bound to apply the dominated convergence theorem. By 
definition of Ojy, Fri[Hi\j > 9^] < e^^ for all 77 £ n^. By the Markov property, we 
obtain that P,,[Hn > t§N] < fr,[HN > M^jv] < e^W so that P^^[Hn/0n > t] < 
e-[*I. D 

Corollary 3.6. Assume that the hypotheses of Theorem \3.1\ are fulfilled. Let 
{fiN '■ N > 1} be a sequence of probability measures and suppose that there ex- 
ists a sequence Sn, T^^^ ^ Sn ^ Ey„[i/jv], such that 

lim P^^[Hn<Sn] = 0. (3.9) 

N—>-oo 

Then, under ^n, i/Ar/E^„[iJjv] converges in distribution to a mean one exponential 
random variable. 

Proof. Let Uat = E^„ [Hm] and fix i > 0. Clearly, 

P^J.AHN<tVN] = P^ASN<HN<tVN] + P^ASN>HN.HN<tVN] ■ 

By assumption, the second term on the right hand side vanishes as TV f 00 1 while 
the first one, by the Markov property, is equal to 

E^„ [ ll^AT < i/^} P^(S„) [ i^AT < t Uat - 5a. ] ] 
= ^^^[PriiS^HN <t\]N - Sn]] 

- ^^^[1{Sn > HN}V-,-iiSr,)[HN <tUN - Sn]] . 

As before, the second term on right hand side vanishes as N ^ 00. The first one, 
since T™"^ ^ Sn is equal to 

IPi^jv [ Hn < t Vn — Sn] + Rn , 

where limjv -Rtv = 0. Since Sn "^ Utv, by Theorem 13.11 the first term in the 
previous displayed formula converges to 1 — e^*, which proves the corollary. D 

To apply the previous corollary one needs among other things to estimate P^„ [ Hn < 
Sn ] and Ejy„ [Hn] ■ In the next section we present a general method to estimate the 
latter sequence when a dynamical large deviations principle is available. There are 
several ways to bound Fi^^[Hn < Sn]- We present below three approaches. The 
first two uses the enlarged processes introduced by Bianchi and Gaudilliere [11] , the 
second and the third ones are taken from the martingale approach to metastability 

Consider a sequence 7jv of positive real numbers. Let il^ be a copy of the set 
V,N and recall from [S] Section 2.C] the definition of the enlarged process associated 
to the sequence 77V, a Markov process, denoted by r]*{t), on Hn U $7^ which jumps 
from a state i] G Hn to its copy rj* G fl^ at rate jn- Denote by z/]^ the stationary 
measure of the enlarged process and recall that v'^irj) = Vn{v*) = (l/2)^'jv(??)- Let 
cap^^ be the capacity with respect to the enlarged process, and for a subset B of 
r^Af, denote by B* the copy of the set B. Next result is CoroUary 4.2 in ^. 
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Lemma 3.7. Let ^n be a sequence of probability measures concentrated on A 



and set jn = Sj^ . A 



N 



that 



lim SnE^^ 






cap,(Ajv, (A^)*) = 



Then, (13.91) holds 



Theorems 2.4 and 2.7 in [24] provide variational formulae for the capacity. The 
second theorem expresses the capacity as an infimum over flows. It permits, in 
particular, to obtain simple upper bounds. An elementary bound for the capacity 
is obtained as follows. By definition of the capacity and since v*{ri) = (1/2)2^(77), 

capMNAA'kr) = (1/2) Y. Hv){RN{v,^N) + lN}P;[H^A'j,y <H+J , 

n£AN 

where P^ represent the distribution of the enlarged process r]*{t) starting from rj. 
Therefore, I^M holds if 



lim Sm E, 



N- 



N i^UN 



V dviq / . 



{s,' 



xa&yi RN{r],VtN)] vn{An) ^ Q ■ (3.10) 



Lemma 3.8. Let fXN be a sequence of probability measures on fl^ and let R' {rj, Aiq) — 
\{v^A%}RN{r],AN). Ass^ 



that 



lim |mw(A^) + E,^[R'^{ri,AN)\ ^ Mw(r?) , „, . | = 

iv^oo L J ^ cap^(77,rj^)J 



for some sequence 7^ ^ Sjq. Then, 



holds. 



In the reversible case, the Thomson principle permits to estimate from below 
the capacity. If 7^^ ^ T^^^, starting from any state in fl^, the distribution of 
r]*{Hn* ) is close to the stationary state i^n lifted to r2J^. Since the capacity can be 
interpreted as the inverse of a distance, the sum on the right hand side measures 
the distance from /ijv to the stationary state vn- 



Proof of Lemma \3.8[ We first replace the deterministic sequence Sn in (|3.9p by a 
sequence of exponential random variables independent of the Markov process ri{t). 
Denote by Cn a mean 7^ exponential time independent of the Markov process 
77(t). Since Sn < 7w^ 

hm sup P^„ [Hn <Sn] < hm inf P^„ [ i^Ar < e^ ] . (3.11) 

N^oc N-^co 

Repeating the steps which led to (|3.4I) . we obtain that 



[Hn <cn] < f^NiAN) + Ep 



R'j^{ri{s),AN)ds 



In this step we used twice the monotone convergence theorem and we replaced zn 
by Cat A i to overcome the unboundedness of Zn- 

Clearly, starting from any configuration in J7jv, we may interpret ejv as the 
hitting time of J7^ for the enlarged process so that 



i?^(r/(s),AAr)ds 



R'j^{r]*{s),AN)ds 



12 



O. BENOIS, C. LANDIM, M. MOURRAGUI 



By [71 Proposition A. 2], since the equilibrium potential is bounded by 1 and since 
'^Niv) = (l/2)i^Ar(?7), V £ rtN, the previous expectation is equal to 



Yl ^^N{v)K 



r;eOj\ 



R'N{'n*{s),AN)ds 






which proves the lemma. 



D 



We conclude this section with a third set of sufficient conditions for (13.91) . Denote 



by T^' the relaxation time, i.e. the inverse of the spectral gap of the symmetric 
part of the generator, and denote by || • \\p the norm of LP{i'n), < p < oo. 

Lemma 3.9. Let Sn be an increasing sequence and let fi^ be a sequence of prob- 
ability measures on flN- Assume that 



lim 

N-¥oo 



N- 



|AtJv(AAr) + SNiyN{A%)rN{A%,AN)^ = 0, 
lim \\R'j,i-,Ar,)h "^^"^ 



di^N 



W I 1 - e 



-Sn/T}^ \ _ 



0, 



where _R^(?7, Ajv) = 1{?7 ^ An} RN{r/, An). Then, p.9p holds. 

We may estimate \\R'j^{- ,An)\\1 by ||i?^( • , Ajv)||oo ||-R^( • '^Af)lli and recall 
that ||i?^( • , An) 111 — ^Ar(A^) rAr(v4^, An) which vanishes asymptotically. 

Proof of Lemma \3.9l Repeating the steps which led to p.4p , we obtain that 



"^^[Hn <S'Ar] < fiN{AN) + 



1(77(5) ^ An} Rn{Ms), An) ds 



mi^N) + SNiyN{A'N)rNiA'N,AN) + 



E„ 



RNivis)) 



ds , 



where RnIv) is the j/Ar-mean zero function RNiv) — R'n^'H^ An) — -E-i/jv [R'n^''^ An)], 
and E^^[R'j^{r],AN)] = VN{A%)rN{A%,AN). 

We estimate the last term of the previous displayed equation. Let fsiv), rj e Hn, 
s > 0, be the unique solution of 



foiv) 



fj-Ni-n) 



"/s — LnJs 



r \ , , JS — ^NJs 1 

VNVn) ds 



where L^ stands for the adjoint of Ln in L'^{vn). With this notation the integral 
in the penultimate displayed equation becomes 

Sn rSN 

{RN,fs)uNds < \\RNh / {fs;fs)l';^ds, 
Jo 

where ( • , • )j/„ represents the scalar product in L'^[i'n) and {fs\ fs)vN the variance 

of/s. It is well known that (/s; /s)i,„ < (fo] fo)i^N^^^'^^'^" ■ The previous expression 

is thus bounded by 

\\RNh{fo;fo)li'n^'{l^e-^^/^'^''), 

which proves the lemma by replacing variances by L^ norms. D 
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4. Expectations of hitting times 

We showed in the previous section that in the context of finite state Markov 
processes, the hitting time of rare events is asymptotically distributed according 
to an exponential law. We show in this section that the expectation under the 
stationary measure of these hitting times can be estimated if one is able to prove a 
dynamical large deviations principle. Instead of presenting this result in a general 
setting, we examine the case of the BDSSEP. 

The dynamical large deviation principle. We recall a result first proved in '5] , 
and then in 10 in the form presented below. We say that sequence of configurations 
{rj^ : N > 1}, T]^ G ^n, is associated to the macroscopic density profile p € M if 
the sequence tt-'^ {ij^) converges to p in M as A^ — ?► oo. 

Given T > 0, we denote by D[[0, T];M) the Skorohod space of paths from [0, T] 
to M equipped with its Borel cr-algebra. Elements of D([0,T],M) will be denoted 
by u(t) and sometimes ut- 

Fix a profile 7 G M and consider a sequence {rj^ ■ N > 1} associated to 7. It has 
been proven in ^7\ [27] following the work of [HI [28] that as iV — > 00 the sequence 
of random variables 

n^(t) := n^irj^itN^)) , (4.1) 

which take values in D([0,T],lVi], converges in probability to the unique weak 
solution u{t) of the heat equation with Dirichlet boundary conditions: 

dtu = (l/2)Au , 

u{t,0)=a, u(i, l) = /3, t>0, (4.2) 

u{0,x) =j{x) , 0<a;<l. 

Note that time has been speeded- up by N'^ in (|4.ip . 

Recall the definition of the rate functional ^[o,t]('|7) of the dynamical large 
deviations principle introduced in (j2.3p . The next two results have been proven in 

m- 

Lemma 4.1. Fix 7 G M and T > 0. The functional /[o.t]('|7) ^s lower semicontin- 
uous and has compact level sets. Any path u with finite rate function, /[o,t](''^|7) < 
00, is continuous in time and satisfies the boundary conditions u(0, •) — ^{■), 
u(-,0) — a, w(-,l) — /3. Furthermore, any trajectory u with finite rate function 
can he approximated by a sequence of smooth trajectories {m" : n > 1} in such a 
way that /[g j-](w"|7) converges to /[o,t](u|7)- 

The dynamical large deviation principle can now be stated. 

Theorem 4.2. Fix T > and an initial profile 7 in M. Consider a sequence {ry : 
-^ ^ 1} of configurations associated to 7. Then, the sequence of probability measures 
{P^ o [tt^ {N"^ ■ )^^ : N > 1} on Z3([0, T], M) satisfies a large deviation principle 

with speed N and good rate function IiQrp]i'\l)- Namely, /[o.t]('|7) • -C([0,r];M) -^ 
[0, 00] has compact level sets and for each closed set £ C -D([0, T], M) and each open 

setO C i:i([o,r],M) 

lini sup ^ log P,^„ {tt^ iN^-)e€) < - M /[o^t] {u \ 7) 



Af- 



lim inf - log P^. (tt^ {N^-)eO) > - inf /[o,t] ("I7) 
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The static large deviation principle. The large deviations principle for the 

Q,/3' 



empirical measure under the stationary state i'^ a, stated below, is taken from 



[HUH]. 

Theorem 4.3. The sequence of probability measures {i^a 0°i'^^)~^ : N > 1} onM. 
satisfies a large deviation principle with speed N and good rate function V. Namely, 
y : M — )■ [0, oo] has compact level sets and for each closed set C C M and each open 
setO cM 



limsup-log<^(7r^ee) < -i£V{j) 
1 



N- 



N " "'^^ ^-J - ^ge 



Expectation of hitting times. The main result of this section can now be stated. 
Fix an open subset of M and let 

An - (TT^'r^O) = {77 G f^AT : ^^(ry) e 0} , 

and let Hpf = H^j^ be the hitting time of the set A^- Note that H^ coincides with 
the hitting time Hq introduced in Theorem 12. II 

Theorem 4.4. Fix an open subset of M. Assume that T™"' ^ expjaiV} for 
all a > 0, that Hn/K^n [Hn] converges in distribution to a mean one exponential 
random variable, and that 

V{0) := inf ^(7) = inf ^(7) . 

Then, for every e > 0, 

E^N [Hn] El,™ [Hn] 



In particular. 



hm 1-logE,. [Hn] = y(0) 



To prove this result we first need a dynamical large deviations principle starting 
from the stationary measure. 

Theorem 4.5. For each T > 0, each closed set € C D{[0,T],JA) and each open 
setD CD([0,T],M), 



1 



limsup-logP^^ {TT^{N^-)e€) < -mi_{l[o,T]{u\uo) + V{uo)} 



liminf 1 logP;:'™^ (7r^(A^2 ■)eD) > - mf^ {/[o^t] ("I«o) + ^("o)} • 

Proof. In order to simplify the expressions, we will use the fact that concerning the 
SSEP process, as mentioned in [TOl last part of section 2], the two dynamical rate 
functionals /[o^t]("|7) and /[o^t](^|7) (see (|2.3p ) are the same. 

We start with the proof of the upper bound. The arguments closely follow 
the ones presented in [10] . Theorem 14.31 is used afterwards to estimate the large 
deviations from the initial stationary distribution. 
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It is well known that using an exponential tightness argument, it is enough 
to prove the upper bound for compact sets. For any function {t,x) n> Ht{x) £ 
Cq' ([0,T] X [0, 1]), we introduce the exponential martingale M^ defined by 



Mf = exp{iv[(^f,i?t}-K^,iJo) 



N 



-JV(.f,ff.)(5^_^^^)gJV(.f,H.)rfs]| 



Using a super-exponential estimate ([ini Theorem 3.2]), for any (5 > and e > 0, 



H,N 



there exists a set of configurations rj € Bg ^ such that for any S > 



lim lim sup -— log P^ 



Kff 



and on which 



Mf = expN{jH{TT^^'\n^) + Onie) + 0{S)} 



where the functional Jh was defined in p.2p . Oni^) (resp. 0(d)) is an deterministic 
expression which vanishes as e J, (resp. 5 J, 0) and where, for any density tt G M, 

1 f 

7r'(u) = — / Tr{u)du. 

2^ J[M-£,u+e]n[0,l] 

Let X be a compact subset of D([0,T], M), then 



limsup^logP;^„ [n^ £%] 



1 



< Hm sup Hm sup Hm sup — log P^„ [{tt^ £ X} n Bff] 



and we can write 



^.11 



[{n^ £ %} n B^f] = E^„ [M^{M^)-\^ 



"e3C}nBf;"J 



Therefore, 

— lop-P^ 

N 



logF^.J{n^£X}nB^f] 
1 



< 



N 



logE^« Mf expA^sup{- Jff(u>^)} +0H{e) + 0{6) 



«G3C 



and since M^^ is a mean 1 martingale, we get 



lim sup — log Pf)'™ [tt^ £X] 



< lim sup lim sup — log i?j,jv expiVsup { — JH(u^|7r^)} 



e^O TV^oo 



N 



ueJC 



We notice that the map tt i— >■ sup„gg(;{— Ji:/(u'^|7r)} is continuous on M, so we can 
apply Varadhan's Lemma to the large deviation principle stated in Theorem [ 



lim — log E,,N 



expiVsupj - Ji/(w"|7r^)} 



mSK 



sup { sup{-Jh(u^|7)} - y(7)} = - inf {JH(u'\-i) + ^(7)}. 
7eM «e3C 7eM,ue3C 



^i?^"^^lv^°sP^iv [tt^ e B[o,T]{u,6)] > -/[_5,o](7r|p) - /[o,t](u|mo) 
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Now, since M x 3C is compact, we can follow step by step the arguments of [lOl 
section 3.3] and we get 

lim sup 1 log P;r« W"" eX]<- inf {/[o,t] {u, 7) + ^(7)}, 

which is precisely the required upper bound since /[o,t] ("7 7) < +00 implies that 

Wo = 7- 

The proof of the lower bound is easier. Indeed recalling the definition of the rate 
function V, we only have to show that for any u G D([0,T],M), any S > 0, any 
n G D([—S, 0], M) such that 7r_5 = p and ttq = uq, and for any (5 > 0, 
1 

where BtQ^p-, {u, 5) is the ball centered at u with radius 5 for the Skorohod topology 
on _D([0,T],!M). If we denote by ii the density path given by tt on [— S", 0] and 
u on [0,T], then w e I?([-S',r],M) and I[_s,t]{u) = I[-s,o\{Ap) + ^[o,t]("|wo)- 
Therefore, since v^ a is a stationary distribution, we have 

^<, ['^'^ e ^[o,T](«, -5)] > ^^.^^ [tt^ e S[-s,T](", ^)] • 

As under v^ o the initial empirical density n^ converges to the stationary density 
p, the lower bound proved in [TU] applies here and we get 

hminfi-logP^M [^^eS[_s.T]K'5)] > ~I[-s,T\{u\p)- 

N —¥00 iV a,/3 

D 

Next lemma is also needed in the proof of Theorem 14.41 

Lemma 4.6. Fix a subset 1) ofM andT > 0. Let 21= {u S C([0,T],M) : u{t) e S 
for some < i < T}. Then, 

ini {l[o^T]iu\uo) + V{uo)} > iniV{p). 

Proof. Fix e > and m e 21. Assume that u{ta) G S, < to < T. By ^^, there 
exists To > and a path v G C([-To,0],M) such that v{-Ta) = p, i;(0) = m(0) = 
Mo, -^[-To,o](v|p) ^ l^(Mo) + e- Defining the path w in C[— To,io],^) by z«(t) = v{t), 
—To < t <0, w{t) = u{t), <t < to, we obtain a path connecting p to u{to) e 23. 
By (|2.4p . /[_rn^tQ](w|p) > infpgs ^^(p). It follows from the estimates just obtained 
that 

I[o,T]{u\uo) + V{uo) > I[o^t„]{u\uo) + V{uo) 

pfc -D 

which proves the lemma. D 

Proof of Theorem \4-.4\ Fix e > 0. There exists 7 G such that 

1/(7) < inf ^^(P) + (e/2), 

and there exists S such that Bsij) C 0. By (|2.4p and by translation invariance 
of the dynamical rate function, there exist Te > and a path zt^'^'(i), < t < Tg, 
-u^"' = p, u(")(rc) = 7 such that 

/[o,T.](w^'^|p) < inf F(p) + e. (4.3) 
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For if > 0, T > and a path u G D{[0,T],M) denote by M^^xiu) the open ball 
in Z3([0,r],M) of radius ip centered around u. Let G = B^^t^C"" ); G^ = {w G 
D{[0,LTJ^'''/N^ +T,],M) : u{LT'^'''/N^ + • ) e G} It is clear from the definition 
of G that Gl C {Hn < LTg''' + T,N^}. Hence, for any configuration ^ £ I^at, 

^i[HN < Ln^^^ + T,N'] > ¥^[7T^ e Gl] = ^ PLTjiMtOM^" ^ G] , 

CeHw 

where Pt{ri,(,), t > 0, stands for the transition probability of the BDSSEP. By 
definition of the mixing time, the previous expression is bounded below by 

P,«^ [tt^ e G] - 2-^ . 

Therefore, for every L > 1, 

inf F^[Hn < LTj^'"" + T,N^] > P^iv [tt^ e G] - 2'^ . (4.4) 

By Theorem 14.51 by definition of G, by (|4.3I) and since u^ = p, V{p) — 0, 
lim^nf 1 logP,«^ [tt^ e G] > - mf^ {/[q.t,] ("I^o) + l^(wo)} 

Hence, there exists iVo = N(){e^ 5) such that for all N > Nq, 

P^«.[7r^eG] > exp-iV{y(0)+2e} . 
The previous estimate together with (|4.4p ioi L = £ N gives that for all N > Nq, 






Iterating this estimate M times, gives by the Markov property that 

r 1 / , , ^ , \ M 

Hn > M{INT^''' + T^N^] < (1 - e 



max 



-A'[y(0)+2e] I 9-fAf 



Taking £ large enough and setting M ~ exp{A^[y(0) + 2e]}, we conclude that 



lim sup max Pj 



Hn > e^[^(0H2^1 {£NT^'^ + T,N^} 



< 1 



Since, by assumption, NT^^^ < expjeA^} for N sufficiently large and since we 
assumed that H^ /^i,n [Hn\ converges to a mean one exponential random variable, 
we have that 

gAr[y(0)+3e] 

lim inf > . 

N^oo Ej,jv [Hn] 

Conversely, for fc > 0, let 21^ = {w G D{[kT,, {k + l)r,],M) : u{t) e for some 
kT^<t<{k + 1)TJ. By definition, for every L>1 

L-l 

P.«^ [Hn <LN^T,] < J2 ^.-, l^'' e 2tfe] < L P,«^ [tt" G 21] , 

fe=0 

where 21 = SIq and 21 stands for the closure of 2t. 
By Theorem 1431 

limsup^logP.^jTT^eSt] < - mi_{l[o,T,]{'u\ua) + V{uo)} . 
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By Lemma |4.H we may restrict the supremum to paths u m C'{[0,Tf],M.). In this 
case, 2t is contained on the closed set 21' = {u G C([0,Te],M) : u{t) e for some 
< t < TJ, so that 

hmsup-logP,« [^^G2l] < - inf {/[o,T,]("|wo)+^K)} • 

By Lemma 14.61 infuG2i' {llo,Tc]{'^\'^o) + ^(""o)} ^ ^^^^neo^ip) ^^'^ this latter quan- 
tity is by assumption equal to inf^go ^(p)- Hence, there exists Nq such that for all 
N > No, 

V^N [tt^ G al < exp -N\ inf V(p) ~ e] . 

Taking L — (1/2) exp A^{F(0) — e} we deduce from the previous estimates that 

for N sufficiently large. Since, by assumption, Hpf/E^^N [i?jv] converges in distri- 
bution to a mean one exponential random variable, we conclude from this inequality 
that 

gJV{y(0)-2e} 

^™^^pir — nrr < °° ■ 

D 

5. Hitting times of rare events in BDSSEP 

We prove in this section Theorem l2.1l Denote by Rn{v^ the rate at which the 
BDSSEP ry(i) jumps from r] to ^. Recall from (j2.ip the distance d introduced in M. 
With this choice, by Schwarz inequality, 

d(7,7') < Il7-7'll2, (5.1) 

where || • ||2 stands for the L2 norm. 

Lemma 5.1. Fix an open subset (!) 0/ M such that d{p,0) > 0. Denote by An 
the set of configurations in 51 jv for which tt^ {r]) belongs to 0: An = {77 G il^v ■ 
TT^ [rj) G 0}. Then, there exists a > such that 

rjv(A^,Ajv) < e^'^^ and vl^iAj,) < e""^^ 
for N sufficiently large. 

Proof Let 0^ is the closed set defined by O5 = {7 G M : d(7,0) <5},5> 0. We 
claim that there exists ^ > such that 

inf Vh) > . (5.2) 

7605 

Indeed, let 25 ~ d{p, 0) > 0. It is clear from the definition of Oa that d{p,^) > 6 
for all 7 G Os- On the other hand, by '9, Theorem A.l], 



V{p) > 



iH'-^i-i'-wi'-^^}- 



Therefore, since < a < p{x) < /3 < 1 and in view of (|5.ip . there exists cq > 
such that for all 7 G 0^, 

V(j) > CO / {7(x)-p(x)}'dx > cod(7,p)2 > coS\ 
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Denote by OAn the outer boundary of An ■ 

N-2 

dAN = \J {^^An: a-'^+^e e An} (J {^ ^ An : a'i £ An] . 

x=l z = l,N-l 

Since X^feOAf ^^{11,0 < N, by definition of the average rate rAf(A^, A^r), 
rN{A'i,,AN) < ^ ] N ,y^,pidAN) . 

It is clear that for each (5 > 0, OAn C {r; G i^N ■ '^'^ iv) ^ ^s} for N large enough. 
Hence, by Theorem 14.31 and by (|5.2p . there exists a > such that 

for N sufficiently large. The same bound holds for An, which proves the first part 
of the lemma. D 

Estimation of the mixing time in the BDSSEP. We show in this subsection 
by a coupling argument that 

T™'^ < {l/2)N^ . (5.3) 

This bound is not sharp but sufficient for our purposes. 

Assume that a coupling (?7t, 6) has been defined in the product space V,n x Qn- 
This means that both coordinates evolve has the original BDSSEP and that the pair 
does not leave the diagonal once it reaches it. We denote by P,,^^ the distribution 
of the coupling when the initial configuration is {r],£,). Denote by Ho the couphng 
time, the time the process reaches the diagonal. It is well known that 

Tn'"" < inf {t : max P,^,^[Ho >t]< 1/4} . 

The coupling of two copies of the BDSSEP is defined as follows. Fix two con- 
figurations ?7, ^ in flN- We assume that the particles evolve according to a stirring 
dynamics and that particles are created simultaneously in both coordinates at the 
boundary. In particular, the coupled process has reached the diagonal when all 
initial particles have left the system. Denote by Hj the time the particle initially 
at j E An leaves the system. If there are no particles at j set Hj = and note that 
if j is occupied by an 77-particle and a ^-particle they both leave the system at the 
same time due to the stirring dynamics. With this notation. Ho < max-,- Hj and 
for all t > 

jeAjv 
Under the stirring dynamics, the particle at j performs a symmetric random walk 
until it reaches the boundary. If we denote by H-^ the hitting time of the boundary, 
it is known that Ej[iJ-|-] = (l/2)j(7V - j) < N'^/8. The previous sum is thus 
bounded by N^ /8t, which proves claim ()5.3p . 

Proof of Theorem \2.1\ The first assertion of the proposition follows from Lemma 
15.11 (|5.3p and Theorem 13. II The second one follows from Theorem 14.41 

To prove the third assertion, let ^n — N^ and consider the enlarged process 
associated to this sequence. By (|5.3[) and by the second assertion of the theorem, 
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Since dfiN/d,^^^ = l{v e Bn}vIp{Bn)-\ E,.^^[{d^N /dv^^^f] = <^(SAr)-i. 
Hence, as RNiVi ^n) < N, the expression appearing on the left hand side of p.lOp 
is bounded by N^v^ JAm)/v^ JB^)- By the static large deviation principle, 

liinsupllogj.^^(AAr) < limsupllogi.^^(^^eO) < - inf^\/(7) , 



-yeo 



1 . M ,^ 1 



liniinf^logz.^^(BAr) > liniinf-logz.^^(^^GS°) > - inf ^(7) . 

Therefore, by assumption (P3]) . N'^v^ p{An)/v^ piB^) vanishes as TV t oo- By 
remark (I3.10p . condition p.9p is fulfilled. By Lemma l5.ll and by (|5.3I) . the as- 
sumptions of Theorem 13.11 are in force. The third assertion of the theorem follows 
therefore from Corollarv 13.61 D 
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